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sually a student learns to solve a system of linear equations in two ways: 



‘substitution’ and ‘elimination’. While the two methods will of course 


lead to the same answer they are considered different because the thinking 
process is different. In this paper we solve a system in these two ways to 
demonstrate the similarity in the computation. We then see that changing 
the point of view leads us to a ‘simpler’ way to solve a system of equations. 
This leads naturally to two other consequences, viz., what is known as Chio’s 
pivotal condensation process for computing determinants and Cramer’s Rule. 
While the condensation process for computing determinants is known, it 
is not widely known, and the manner of solving equations developed here 
has not been seen elsewhere. This should be of interest to anyone teaching 
solving systems of linear equations (especially by hand) and can be the basis 
for teaching the basics of solving systems of equations, or for use as a guided 
project. This material is particularly relevant for the topic of matrices in unit 
2 of the Specialist Mathematics, the topic of Algebra and matrices in unit 1 of 
the General Mathematics curriculum, as well as anywhere where multivariate 
applications appear such as finding regression lines in the data collection 
topic in unit 3 of essential mathematics (perhaps as a special project) and the 
bivariate data analysis topic of unit 3 of the general mathematics curriculum. 

Solving systems of linear equations 

While we can deal with any size system of equations, we will begin by considering 
the system of equations given by: 


2x + 3y + z = 4 


< x + 2y-z = -l 
-x + 2y + 2z = 3 

We can solve this in two ways, viz., ‘substitution’ and ‘elimination’. Let us 
consider one at a time. 
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Substitution 

Let us solve each of the equations for x. We get the system 

—3) —z + 4 

x =- 

2 

- x = -2y + z —1 

-2)-2z + 3 

x = —-- 

-1 

Substituting the right-hand side (RHS) of the first equation in for in the 
second and third equation and cross multiplying we find: 

2x + 3y + z = 4 

■ 2(-2y + z-l) = l(-3y-z + 4) 

2(-2) - 2z + 3) = -l(-3y - z + 4) 

or, by multiplying both sides of the second and third equation by (-1) 
(distributing through the parentheses) and collecting variables on the left- 
hand side, we have 

2x + 3y + z = 4 

■ (2 ■ 2 -1 • 3)y + ((2 ■ (-1)) -1 • l)z = 2 ■ (-1) -1-4 

(2-2-(-l)-3)y + (2-2-(-l)-l)z = 2-3-(-l)-4 


or, by recognizing the coefficients as determinants of 2 x 2 matrices, we have 



2x + 3) + z 

= 4 


2 3 

\\ | 

2 1 


2 4 

1 2 

) + 

1 -1 

z — 

1 -1 

2 3 

1 

2 1 


2 4 

-1 2 

y + 

-1 2 

z — 

-1 3 


or, 


2x + 3y + z = 4 
< y - 3z = -6 
7y + 5z = 10 

So now each of the second and third equation can be solved for and cross 
multiplying yields: 


2x + 3y + z = 4 
y = 3z - 6 
-5z + 10 
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or 


2x + Sy + z = 4 
y = 3z-6 

7(3z - 6) = 1 • (-5z + 10) 


So we have 

2x + 3y + z = 4 
- y - 3z = -6 

(1 • 5 - 7 • (-3))z = 110- 7(-6) 


or 


or, 


1 

7 


2x + 3y + z = 4 
y 3z = —6 



1 -6 
7 10 


2x + 3y + z = 4 
< y - 3z = -6 
26z = 52 


This can then be solved (via back substitution) to yield the solution z = 2, y = 0, 
x = 1. 


There is another option for this problem which we now describe. Going back 
to 

2x + 3y + z = 4 
y = 3z - 6 

< 

-5z + 10 
. ) ' = ^ — 

we see that we can extend this substitution step to include the first equation 
(substituting the value of y found from the second equation in both the first 
and the third): 

—2x - z + 4 

y = - 

3 

< y = 3z-6 
-5z + 10 
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to get 

3 (3z - 6) = 1 • (~2x - z + 4) 

< y = 3z - 6 

7 (3z - 6) = 1 ■ (-5z + 10) 

or, by multiplying both sides of the first and third equation by -1 (distributing 
that -1 through the parentheses) and collecting variables, we have 
(2 • 1 - 0 • 3)x + (l ■ 1 - 3 • (-3)) z = 1 • 4 - 3 • (-6) 

< y = 3z-6 

(l ■ 5 - 7 • (-3))z = 1 ■ 10 - 7 • (-6) 


or, by recognizing the coefficients as determinants of 2 x 2 matrices, we have 


2 3 
0 1 


x — 


3 1 
1 -3 


z = — 


3 4 

1 -6 


y — 3z = -6 


1 -3 
7 5 


z 


1 -6 
7 10 


or, 


2x + 10z = 22 
y — 3z = -6 
26z = 52 


Now substituting in for z in both the first and second equation: 

-2x + 22 


z = ■ 


So, 


10 
= z y-6 
-3 
_ 52 
Z _ 26 


52 ■ 10 = 26(-2x + 22) 
52(-3) = 26(—> - 6) 
52 


z = - 


26 
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or, 


2 

10 

X = — 

10 

22 

0 

26 


26 

52 

1 

-3 


-3 

-6 

0 

26 

y=- 

26 

52 


26z = 52 


52x = 52 
< 26 y = 0 

26z = 52 


from which we can readily see that, as we found earlier, x = 1, y = 0 and z = 2. 

Elimination 

Now let us solve the same problem using elimination (boxing what is known 
as the pivot). We begin with the equation 


' ® 

3 

1 

\ 

4 

1 

2 

-1 

-1 

l - 1 

2 

2 

3 J 


and see that we should replace the second row (R2) by twice row 2 minus row 
1 (2 x R2 - Rl) and the third row (R3) by twice row three minus negative one 
times row one (2 X R3 - (-l)Rl). Of course, 2 x R3 - (-l)Rl = 2 x R3 + Rl, but 
the formulation we use is more algorithmic in that it allows for easy adaptation 
to a different problem by replacing 2 and -1 by appropriate numbers taken 
from the matrix for that problem. So the new row two is 

( (21-1-2) (2-2 —1-3) (2•(-!)-1-1) (2-(-l) —1-4) ) 


f 

0 

V 

2 

3 


2 

1 


2 

4 

\ 

1 

2 


1 

-1 


1 

3 

) 


and the new row three is 

( (2•(-!)•-(-!)• 2) (2-2-(-l)-3) (2-2-(-l)-l) (2-3-(-l)-4) ) 


f 

0 

V 

2 

3 


2 

1 


2 

4 

\ 

-1 

2 


-1 

2 


-1 

3 

) 


47 


On solving systems of equations by successive reduction using 2*2 matrices Australian Senior Mathematics Journal vol. 28 no. 1 
































Australian Senior Mathematics Journal vol. 28 no. 1 Carley 


So we have 


f 1—1 



A 


2 

3 

® 

3 

1 

4 



2 3 

1 

2 

-1 

-1 

= 

0 

1 2 

l - 1 

2 

2 

3 

7 




2 3 

-1 2 


2 3 1 

0 0 -3 
0 7 5 


1 

2 1 
1 -1 

2 1 
-1 2 


4 

2 4 
1 -1 

2 4 
-1 3 


4 

-6 

10 


Now using the boxed 1 as the pivot we want to replace R3 by (1 x R3 - 7 x R2) 
and (optionally if back-substitution will be used) R1 by (1 x R1 - 3 x R2), i.e., 
the new R3 are R1 are 


Rl 


= ( (1-2-0-3) (1-3-3 1) (1 • 1 -(-3)• 3) (1 • 4-(-6)• 3) ) 


2 3 
0 1 


3 1 
1 -3 


3 4 

1 -6 


and 


i?3 = ( (10-0-7) (1-7-7-1) (1-5-7-(-3)) (l-10-7-(-6)) ) 


0 1 
0 7 


1 -3 
7 5 


1 -6 
7 10 


giving 


' 2 

3 

1 

4 

0 

1 

-3 

-6 

l 0 

7 

5 

10 J 


2 

3 

1 

o 

3 1 

1 

1 

1 

3 4 

0 

1 


1 -3 

1 

1 

1 

1 -6 


0 

0 


1 

0 


-3 

1 -3 
7 5 


-6 

1 -6 
7 10 


2 

0 

10 

! 22 

0 

1 

-3 

-6 

0 

0 

26 

! 52 


While at this point we could solve using back substitution, we can also use 
the boxed value as the pivot we can replace Rl with (26 x Rl - 10 x R3) and 
R2 with (26 x R2 - (-3) x R3) giving the new rows: 

f?l = (- (26-2-0-10) 0 (26-10-10-26) (26-22-10-52) ) 


2 10 
0 26 


0 10 
0 26 


0 - 


10 22 
26 52 
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f?2 = (- (26 • 0 - (-3) • 0) —(26 -1 - (-3) - 0) 0 -(26-(-6)-(-3)-52) ) 


f 

0 

-3 


1 

-3 

0 - 

1 

CO 

1 

\ 


0 

26 


0 

26 

26 52 


V 








y 


giving 


' 2 

0 

10 

i 22 "l 

0 

1 

-3 

-6 

0 

0 

26 

; 52 


v — y 


2 10 
0 26 


0 


1 -3 
0 


52 0 0 

0 26 0 
0 0 26 


0 


52 

0 

52 



10 22 

A 


26 52 



-3 -6 



26 52 



52 

y 


from which we can again see that x = 1, y = 0 and z = 2. 


Comments on the two methods and a change of the 
point of view 


Both of the above methods yield the following calculation in matrix notation: 


GE<] 

3 1 

A 

4 


2 

3 

2 3 

l 

l - 1 

2 -1 

2 2 

-1 

3 

2 


0 

1 2 


2 

-1 


1 

2 1 
1 -1 

2 1 
-1 2 


2 4 
1 -1 

2 4 
-1 3 


which is equal to 


' 2 

3 

1 

\ 

4 

0 

a 

-3 

-6 

V 0 

7 

5 

10 

y 


0 

0 

2 

0 


3 

1 

0 


1 

-3 


1 -3 

1 

1 

1 

1 -6 

7 5 

i 

1 

1 

7 10 


4 

-6 




y 


3 1 
1 -3 
0 26 


4 

-6 

52 


\ 


y 
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or 


' 2 

3 

1 

\ 

4 

0 

s 

-3 

-6 

. 0 

7 

5 

10 


2 3 

1 

o 

3 

1 

1 

1 

1 

3 4 

0 1 


1 

-3 

1 

1 

1 

1 -6 


-3 


2 0 10 
0 1 -3 
0 0 26 


22 

-6 

52 


1 -3 

I 

1 

1 

1 -6 

7 5 

1 

1 

1 

7 10 


-6 


which can be written as 


' 2 

0 

10 

i 22 "l 

0 

1 

-3 

-6 

0 

0 

26 

; 52 


v — y 


2 10 
0 26 


0 


0 


52 0 

0 26 


1 -3 
0 26 

0 


0 

0 

26 


52 

0 

52 



10 22 

\ 


26 52 



-3 -6 



26 52 



52 

) 


from which we can again see that x = 1, y = 0 and z = 2. 


The points of view of the methods of substitution and elimination are 
different (that is to say that thought processes are different). The substitution 
point of view involves a fair amount of algebraic manipulations to proceed 
from one step to another (and the augmented matrix formulation is unnatural 
using this approach), the elimination point of view involves dealing with adding 
multiples of vectors to each other in order to form the next augmented matrix 
(which is often done on the side). But in both cases the outcome of each step 
can be written simply by replacing certain elements of the augmented matrix 
with 2x2 determinants which are formed by using ‘submatrices’ and minus 
signs according to where the pivot is in relation to the opposite corner of 
the submatrix (the ‘counter pivot’). The sign situation is summarised in the 
following diagram where the pivot is in the centre and the plus or minus signs 
are the locations of the counter pivot: 

( + * - ) 

* pivot * 

- * + 

It should also be noted that there is some freedom in choosing pivots 
(which variables to substitute or which columns to ‘wipe out’) and they should 
be chosen in order that they not be zero (for then there is a division by zero 


50 










































or a non-invertible row replacement, in the substitution and elimination views, 
respectively). 


Example 1 

Here is another example where different pivots are used. 


0 

-1 

1 

\ 

1 


' -2 

0 

0 

-2 " 

2 

1 

ES 

1 

= 

2 

1 

-1 

1 

l 2 

0 

-2 

-2 

2 


2 

m 

0 

4 

) 


ED 0 0 
2 0-2 
2 2 0 

-4 0 0 

0 0 8 
0-8 0 


-4 

-2 

4 


-4 

16 

-8 


\ 


which leads to x = 1, y = 1, and z = 2. 

Note: This procedure works for any size system of linear equations and can 
be used to solve equations of the form Ax = b for various b simultaneously; 
therefore it can easily be used to compute the inverse of a matrix (by starting 
with the identity on the right side of the dotted line). 


Determinants 

Note that a determinant can be computed by using row operations, that is, by 
elimination. This procedure is based on the transformations described above 
and the following properties of determinants: 

. a-\A\ = |A'| (or | A\ = i|A'|), where the ith row has been replaced by a times 
the till row plus b times the /th row. 

a-|A'| 

if A' is obtained from A by exchanging two adjacent rows or 


0 A' 

• I4=-|A1 

columns. 


So the determinant of the matrix 

' 2 

1 

l - ] 

can be computed by using the exact same procedure as above with the 
adjustments coming from the three changes above of which, in this case, only 
the first two are used (the third point is used when other pivots are chosen as 
we will demonstrate with examples below). We perform the calculation here: 


3 

2 

2 


1 

-1 

2 
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( r —I 5 

@ 3 1 

1 


3 

2 3 

det 

1 2 -1 
l - 1 2 2 J 

= —5-det 
2 2 

0 

1 2 


= —n-det 


2 3 

-1 2 


2 3 1 

0 0 -3 
0 7 5 


1 

2 1 

1 -1 

2 1 
-1 2 


Now we will use the 1 as a pivot, but since the determinant of a triangular 
matrix is as easy to compute as that of a diagonal matrix, we will only apply the 
operation with the pivot to row 3: 


-^rdet 


2 3 1 

0 0 -3 
0 7 5 


1 1 , 
^ det 


2 3 

0 [l] 

0 0 


1 

-3 

1 -3 
7 5 


1 1 , 

= ^r^r det 


2 3 1 

0 1 -3 
0 0 26 


= X \-26 = 13 
2-1 1 

Generally, the determinant of an n x n matrix A = ( ajj with a n ^ 0 can be 
written as the reciprocal of a u n ~~ times the determinant of an (n — 1) x (n — 1) 
matrix thusly: 

det(A) = —det(A') 


a 


n-2 

11 


where A' is an (n — 1) x (n — 1) matrix whose ij th entry is 

f ^ 


det 


*u 


i.j+i 


a i+l,l a i+\,i+\ 

V 1 J 


The pivot (which in the above case is a n ) may also be taken to be in the 
/^-position, provided that number is not zero, in which case the above 
formulation is that 

/ -i \ row+column of pivot 

det(A) = ^2 -_-det(A') 


where A' is an (n - 1) x (n - 1) matrix whose entries are \a kl a l} - a h] a il \, 
i ^ k and j l, arranged in the same order as the subscripts ij appear in the 
matrix (since i^k andj i 1 1 there is one less row and one less column in A' than 
in A). Note that a kl a ij - a hj a u is plus or minus (depending on the relationship 
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between kl and ij) the determinant of a ‘sub-matrix’. This method is known 
as Chio’s pivotal condensation method, though according to Howard Eves, 
there are traces of the method in an earlier paper by Hermite (see Eves, 1966, 
p. 129). It is likely that it has been discovered a number of times. 


We give a few more examples to illustrate the procedure (using different 
pivots which are boxed). 

Example 2 

This example demonstrates why a (-1) shows up when using some pivot values 
(note the column swap needed to put the matrix in diagonal form). 

3 2 1 

0 \5\ -5 
0 7 5 

where we have swapped the first two rows to avoid a non-zero pivot. Continuing, 


2 

IS 

1 

1 

2 

3 

1 


3 

2 

-1 

5 

0 

-5 

1 

-2 

1 


7 

0 

5 



= (_ 1 ). J_J_. 

1 > 3 2 5 1 


3 2 1 

0 5-5 
0 0 60 


=< “ 1) W' 60=H) '?'5 r60= “ 20 


Example 3 

This is an example using a larger matrix. 


det 


' IS 

1 

0 

1 

\ 

1 


( 1—1 



\ 

1 

1 

1 

1 

1 

1 , 

= —5- det 

0 

2 

1 

1 

2 

0 

2 

0 

2 

-2 

4 

-2 

2 

0 

1 

3 

3 

0 

2 3 

2 

6 

6 

0 

l 0 

0 

1 

-1 

1 J 


l 0 

2 

-2 

2 J 


1 1 ^ 

= 7^ det 


( 


V 


\8\ 0 4 

2 4-2 

2-2 2 


1 1 1 , 

= —s-^-det 

2 3 l 2 8 1 


32 

-16 


-24 

8 


= 444 s 2 (- 2 ) 

2 3 l 2 8 1 

= -2 
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Example 4 

This example uses various choice of pivots. 


1 

2 

1 

2 

1 









-1 

-1 

1 

-2 

1 

-2 








-1 

-4 

2 

1 

0 

0 

-1 



0 

1 

-1 

1 

1 

(-2) 3 

-1 

1 

1 

1 

1 

1 


-1 


-5 -5 0 

-2 0 [|] 
-1 -1 -4 
-3 -3 0 


-1 


(_ 2 ) 3 2 2 

-1 (- 1) 1+1 


-2 -10 -10 


-14 -10 -2 

-2 -6 -6 


(- 2 ) 3 2 2 (- 2) 3-2 

- A -8(-16) 

(- 2 ) 3 2 2 (- 2 ) 


-120 -136 
-8 -8 


= 2 


One variation of the method includes the application of the property of 
determinants: |A| = a\A'\ if A' is obtained by A by dividing every member of 
one row or column by a. For example, during the previous computation we 
could have proceeded from the second line as follows: 


Example 4b 


-1 1 

Hi) 3 ? 


H 

-14 

-2 


-10 -10 

-10 -2 

-6 -6 


-1 

(— 2 ) 3 2 2 


H) 3 


a 5 

7 5 
1 3 


5 

1 

3 


-1(-1) 1+1 -30 -34 
~ 2 2 l 3 " 2 -2 -2 

= 2 


This method of computing determinants by hand has its advantages. As 
an illustration of one advantage, it should be noted however that in the 3x3 
determinant above, there a multiplier that can not be ‘gotten rid of’ in the 
same way as the (-2) 3 . Similarly, in Example 2, the multiplier of 1/3 could 
not be cancelled by using the properties of determinants. So the fact that our 
final answer is an integer is encouraging. If we had made a mistake along the 
way, we may have arrived at a non-integer answer and been forced to admit 
our mistake! 
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Cramer’s rule 


Let us now discover Cramer’s rule in the context of an example. Consider the 
solution to the same equation as before: 


' ® 

3 

1 

\ 

4 

1 

2 

-1 

-1 

l - 1 

2 

2 

3 J 


We can solve this for z by putting in triangular form: 


/ r-i ' 5 

§ 3 1 | 4 


/ i \ 

2 3 1 | 4 

1 


/ 1 \ 
2 3 114 

1 2 -1 i -1 

= 

0 0 -3 | -6 

= 

0 1 -3 -6 

1-12 2 | 3 J 


07 5 i 10 

V i / 


0 0 26 i 52 

V i J 


So we find that z = 52/26. But the operations on the third and fourth column 
were identical and the same as the operations used to find the determinant of 
the matrix as described in the previous section. That is to say that 


2 3 1 
1 2 -1 
-12 2 


= K- 21-26 


where ft is a number depending on the pivots and their locations. But 


2 3 4 

1 2 -1 
-12 3 


= ff • 2 T • 52 


for the same K. Therefore 




2 

3 

4 



1 

2 

-1 

52 _ 

K 2T-52_ 

-1 

2 

3 

26 

K- 21-26 - 

2 

3 

1 



1 

2 

-1 



-1 

2 

2 


which is what Cramer’s rule says. 


Rearranging columns we see that 


' E 

3 

2 

4 


-i 

2 

1 

-1 

= 

l 2 

2 

-1 

3 J 

V 


So, by the same argument 


3 2 

4 



' 1 

3 

2 i 

4 ' 

0 3 

3 


= 

0 


3 

3 

1 

1 

-5 

J 


, 0 

0 

-13 i 

i 

- 13 , 
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13 _ A ' l - 5 (- 13 ) _ 

1 3 4 

-1 2 -1 

2 2 3 


4 3 1 

-1 2 -1 

3 2 3 

13 A ' l - 5 (- 13 ) _ 

1 3 2 


2 3 1 


-1 2 1 


1 2 -1 


2 2-1 


-12 3 


And finally: 



1 

3 

4 



' 2 

1 

3 

4 

A 


f 2 

1 

3 i 

4 1 

1 

-1 

2 

-1 


= 

0 

-3 

2 

-6 


= 

0 

-3 

2 

-6 

l - 1 

2 

2 

3 



l 0 

5 

7 

10 



l 0 

0 

- 3 ! 

0 J 


So, 


0 _ K "-2 (- 3 ) 0 

2 1 4 

1 -1 -1 

-12 3 


2 4 1 

1 -1 -1 

-13 2 

-31 ~~ K " ■ 2 • (- 3 ) • (- 31 ) _ 

2 1 3 


2 3 1 


1 -1 2 


1 2 -1 


-12 2 


-12 2 


Conclusion 

We have shown in die context of an example that the methods of substitution 
and elimination for solving systems of linear equations lead to the same 
computation, though through different thought processes. If we adjust the 
point of view still again to a component-wise computation, the calculations 
become simpler (no side-work is needed). This directly leads to the discovery 
of what is known as Chio’s pivotal condensation process for computing 
determinants and Cramer’s rule. 
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